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We extend the investigation of the gravitational collapse of a spherically symmetric Yang-Mills 
| field in Einstein gravity and show that, within the black hole regime, a new kind of critical behavior 

Q\ . arises which separates black holes formed via Type I collapse from black holes formed through 

^\ ' Type II collapse. Further, we provide evidence that these new attracting critical solutions are in 

fact the previously discovered colored black holes with a single unstable mode. 
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I. INTRODUCTION 



■ It is by now relatively well known that gravitational collapse can produce rich structure even within highly simplified 
systems such as spherical symmetry. In particular, near the threshold of black hole formation, the strong field dynamics 
of general relativity exhibits critical phenomena. 

The pioneering work demonstrating this critical behavior in the collapse of a single massless scalar field Q has been 
supplemented by investigations of gravitational waves [§, a perfect fluid || and a variety of other matter models, all 
of which exhibit the same general characteristics. Indeed, to our knowledge no system which has been studied in this 
\ context has been shown not to exhibit this critical phenomena. 

At this point in the subject's development, dynamical evolutions (solutions of the full partial differential equations 
^\ ' of motion), in tandem with analytic and perturbative calculations have given us a reasonable understanding of many 
0^ ■ of the phenomenological details of critical behavior in collapse. (See || for an excellent review of the subject.) 
""q ' In light of this, one of the more interesting discoveries in some of the more recently studied models [||-^| is the 
q-i presence of two distinct types of behavior at the threshold of black hole formation. Specifically, in these models, 
certain regions of parameter space (initial-data space) are found to yield near-critical collapsing configurations which 
display self-similarity and, in the super-critical regime, a scaling law which is continuous in the black hole mass. By 
analogy with the theory of phase transitions, this is called a Type II transition. However, it is found that other regions 
of parameter space lead to critical collapse which has a static (or periodic) solution as an intermediate attractor — this 
results in a black hole transition with a nonzero mass gap. Again, in analogy with the discontinuous behavior in order 
parameter which frequently accompanies first order phase transitions, this is called Type I behavior. For the case of the 
Einstein- Yang-Mills (EYM) model studied in ||), the static solution appearing in Type I collapse is the well-known 
n = 1 Bartnik-Mckinnon solution [fTo] , while for the massive scalar model considered in |^|, the Type I threshold 
solutions are apparently unstable members of the family of "oscillating soliton stars" which have previously been 
constructed by Seidel and Suen albeit in a different context. Finally, in the case of Einstein-Skyrme (ES) collapse 
considered in a static Type I solution was observed, which, as in the EYM case, had previously been constructed 
and studied within a purely static ansatz p2[ . As noted above, each of these three models also exhibits Type II 
behavior — in the case of massive scalar collapse || , the Type II critical solution is the same one originally observed 
in massless scalar collapse |5J], and, interestingly, in the ES model, the Type II solution is evidently identical |)| to 
that observed in the EYM model Heuristically, one expects Type II behavior in any collapse model where the 
initial configuration can be made sufficiently "kinetic-energy-dominated" (ultra-relativistic), whereas Type I behavior 
is expected only in those models which have an intrinsic length scale (or equivalently, mass scale), and which have 
some type of self-interaction which can "balance" the attractive gravitational interaction. 



*Electronic address: matt@einstein.ph.utexas.edu 
^Electronic address: ehirsch@bach.liu.edu 
"''Electronic address: marsa@einstein.ph.utexas.edu 



f 



As noted in the concluding remarks of || , hints of further interesting phenomenology in the EYM model have been 
seen in the super-critical regime where all evolutions are characterized by black hole formation. In this paper we 
study this regime in more detail and present evidence for a new type of critical transition in which the intermediate 
attractors are the "minimally unstable" (one unstable mode in perturbation theory) colored black holes discovered 
by Bizon [O], and independently by Volkov and Gal'tsov [Q. This result is, of course, analogous to the discovery 
that the n = 1 Bartnik-Mckinnon solution is the intermediate attractor for Type I collapse. Crucially, in order to 
accurately model the dynamics of super-critical solutions for long times after the formation of an event horizon, we 
use black-hole excision techniques. Such methods were first successfully employed in a dynamical context by Seidel 
and Suen [[l5|, and have subsequently been studied and implemented by many other authors (see and references 
therein for a more extensive discussion). However, to our knowledge, this is the first time that excision has been used 
to study critical collapse, and we feel our results highlight the power and potential of the strategy to elucidate issues 
relating to the formation and long-time evolution of black holes. Our adoption of excising techniques necessitates the 
use of a different coordinate system than that used in — that work used polar slicing and areal spatial coordinates, 
a system which generalizes the usual Schwarzschild coordinates to time-dependent spacetimes. As is well-known, the 
t = constant slices in the polar/areal system cannot penetrate apparent horizons — thus, for all practical purposes, the 
slices remain outside of event horizons and therefore cannot be used in conjunction with excision. We therefore retain 
areal spatial coordinates, but adopt maximal slicing — in this case the slices do cross apparent and event horizons, and 
excision techniques can be used. 

The outline of the remainder of the paper is as follows. The next section describes the EYM model and the equations 
of motion we subsequently solve numerically. We pay particular attention to gauge and coordinate choices, as well as 
to regularity, boundary and initial conditions for the fields. In Section III we describe our numerical scheme, focusing 
on our specific coordinate choices and on some details of our black-hole excising technique. In Sections IV and ^ we 
describe our results and conclusions. 



II. EQUATIONS AND ASSUMPTIONS 

We are interested in investigating the gravitational collapse of a self-gravitating Yang-Mills gauge field. To begin, 
let us consider the action for an EYM theory 
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where F" is the Yang-Mills field strength tensor. On varying the action with respect to the metric, g^, and the 
gauge connection, A", we get the general equations of motion. We simplify these further by making some additional 
assumptions. In particular, we choose the gauge group to be SU (2) and focus on spherically symmetric gravitational 
collapse. This places restrictions on both the spacetime metric and the form of the ga uge connection. Even so, the 
equations we derive have a rather general form as can be seen in Equations (A6-A22). Subsets of these equations 
have been considered in a variety of different contexts. In particular [[| evolved a version of the Einstein- Yang-Mills 
equations in polar (K e g = 0), areal (b — 1) coordinates with the additional assumption on the Yang-Mills field that 
the connection A^ was purely magnetic. 

Since our interest here is to consider the same model studied in || , but to penetrate into the super-critical regime of 
the phase space, we will likewise make the "magnetic ansatz." As described in the Appendix, together with appropriate 
gauge choices, this ansatz effectively sets all but one of the components of the gauge connection to zero. 

In addition to making these gauge choices we must also choose a coordinate system. As mentioned in the Intro- 
duction, in order to evolve the system for long times to the future of black hole formation, we choose maximal time 
slices and areal (or radial) spatial coordinates. 

As detailed in the Appendix, the equations can now be written in the following form. The evolution equations are 



n = 


/fflH 


- _$ 






a 








$ = 








. a 




w = — n + 


0w' 



aa , . 

+ ~2 W ( 1 ~ W ) 



(2) 

(3) 
(4) 



and the constraint equations are 
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where the matter stress-energy terms are given by 
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We also note that we have an algebraic relation for the sole component, j3, of the shift vector, (3 l = (J3, 0, 0): 

= arK e 6 . (11) 

In addition to the equations of motion, we need boundary conditions on the fields. These are determined by 
demanding regularity at the origin of spherical symmetry, and by enforcing an outgoing condition on the radiation 
fields at large radius. This latter condition assumes that there is no radiation coming in from outside our finite 
mesh. This is not completely true, as in general there will be backscattering of the propagating fields off regions 
of high curvature. However, if our domain of integration is large enough, the contributions from this scattering are 
dynamically negligible. 

The demand for a regular origin results in boundary conditions on the fields at r — 0. From equations (^|) and (|?]) 
it can be seen that a'(t, 0) = a'(t, 0) = 0, as well as a(t,0) = 1. In addition, (3(t,0) — K$ e (t,0) = 0. For the matter 
fields, w(t, 0) = ±1 and w'(t,0) = 0, so that the auxiliary variables II and <I> are both zero at the origin. We note 
that these conditions constrain the Yang-Mills field, w, to be in a vaccum state at r — 0. 

For initial data, we choose a time-symmetric kink for the gauge potential which was previously used in ||. This 
pulse is given by 



w(0,r) = 



tanh (12) 



w(0,r) = 0. (13) 

with the parameters a and b chosen such that u>(0,0) = 1 and u/(0,0) = 0. The two parameters x and s define the 
center and width of the pulse respectively. They also serve as the two parameters which we will vary in order to 
explore the phase space. We note that the implementation of this data in (El was incorrect but ultimately had no 
effect on the overall conclusions of that work. We have fixed the implementation of the kink data here and note a 
minor improvement in the conservation of energy. 



III. NUMERICAL APPROACH AND BLACK HOLE EXCISION 



It should be emphasized that we do not incorporate any form of adaptive mesh refinement into our numerical 
approach. Instead, we use a fixed uniform grid with a mesh spacing which is sufficiently fine to uncover the new 
critical behavior. The drawback to this approach, of course, is that we are unable to fully resolve the discretely self 
similar solutions which arise near the Type II black hole threshold. However, our primary interest here is in the 
supercritical regime, and we are satisfied that previous work has established the nature of the Type II transition. 

What is important in this work is the use of black hole excision techniques which allow us to evolve well beyond 
the formation of the black hole. We discretize the evolution equations (|]-f§) using two time levels with centered time 
differences, and angled spatial differences as described in p6j . The constraint equations are then integrated outward 
from the origin. Equation (^|), the slicing equation, is rewritten in first order form 

5> = S ( a - 2 ~) + ^ (a 2 1 + 2 ^L ) + AnGa (S - 8p) (14) 
a' = 5. (15) 
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The apparent horizon equation in maximal-areal coordinates is simply 



arK% = 1. (16) 

Thus, we can compute a characteristic function C such that 

= f 1 if arK 9 g < 
\ if arK e g > fi H . 

where \in is a threshold value, ostensibly unity, but set to a slightly larger value in practice — typically, fin = 1-1 — to 
ensure that the inner boundary lies strictly within the apparent horizon. At any time, the characteristic function tells 
us which grid points are deemed to be inside (C(r, t) — 0) or outside (C(r, t) — 1) the horizon. As the system evolves, 
we recompute C at each time step and simply "discard" any grid points for which C — 0. The radial locations of these 
grid points will then generally satisfy r < r#(t), where rn(t) is the instantaneous inner boundary of the computational 
domain, and roughly coincides with the location of the apparent horizon. We need no boundary conditions for the 
functions evolved with equations since the characteristic directions at r = rjj are such that only events with 

r > ru are in the past domain of dependence. We can therefore use the discrete form of the evolution equations up 
to and including the inner boundary point. However, the constraint equations (J^-^j must have an inner boundary 
condition. Until a horizon is formed, we can simply use conditions derived from regularity. Once a horizon forms, 
however, we need an alternative. If the horizon initially forms at time t — T and radial coordinate r = R, then for 
a function / £ {a, K e g}, we use f(R,T) as the initial boundary value for /. We then evolve / along this boundary 
using the evolution equation for /. The evolution equations that we use are as follows: 
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Since the lapse a does not have an evolution equation, we simply leave a and 8 "frozen" at whatever specific values 
they have when excision at that particular radius, rjj(t), begins. 

Finally, we note that the programs which we used to generate the results that we will describe below were written 
in RNPL (Rapid Numerical Prototyping Language) pjj. This language has been specifically designed to facilitate 
the differencing and subsequent numerical solution of partial differential equations. 



IV. RESULTS 

On evolution of the equations, we verify previously examined aspects of this system as discussed in ||. In particular, 
we find that in different regions of the parameter space, there are both types of critical behavior at the threshold of 
black hole formation. We observe Type I behavior, wherein black holes turn on at finite mass and the intermediate 
attractor is the n = 1 Bartnik-Mckinnon solution. We also confirm aspects of Type II collapse in the appropriate 
region of parameter space. This includes such things as echoing and the scaling relation for the black hole mass. 
However, because our code does not employ mesh-refinement, we are unable to fully reproduce the results of 0. 
However, as discussed previously, our main emphasis here is to address the outstanding question of what happens in 
the super-critical region and to the future of black hole formation. 

On evolving into the supercritical (black hole) regime, we find that the dynamics of the resulting spacctimcs 
exhibit two distinct types of behavior which we describe below. It turns out that these can be characterized by the 
asymptotic vacuum state of the Yang-Mills potential w. These two behaviors are related to the Type I and Type II 
critical collapses identified previously. 

The original Type I collapse is described, in part, by a black hole which forms with finite mass. Slightly super- 
critical Type I evolutions, now continued well to the future of black hole formation confirm this behavior, as one would 
expect. In addition, after the black hole is formed, part of the remaining Yang-Mills field outside the horizon radiates 
to infinity and some of it falls down the black hole causing the black hole to grow. How much larger the initial black 
hole becomes through the infall of additional matter depends, in part, on where in the parameter space the evolution 
begins. Eventually, however, and in accordance with what might be expected from the no-hair theorems, the exterior 
spacetime settles down to a Schwarzschild black hole with the Yang-Mills field in its vacuum state. In general, for 
slightly super-critical collapse (either Type I or Type II), as the remaining mass-energy radiates away the gauge 
potential approaches one of its two vacuum values: w — ±1. For the case of slightly super-critical Type I collapse, 
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in the asymptotic regime, the gauge potential always picks the w = — 1 vacuum state (recall that our initial data 
family (12) has w(t,0) = +1 and w(t, oo) = —1). This turns out to be true even for collapsing configurations which 
are no longer just slightly super-critical but are, in fact, well into the black hole regime. We can thus characterize a 
sort of "generalized" super-critical Type I collapse by the asymptotic value of the gauge potential w. 

When Type II collapse is considered, the overall picture is broadly the same. Again, after the original black hole 
forms (at finite mass since we are slightly super-critical) some of the matter remaining in the region exterior to the 
black hole falls into the black hole causing it to grow and the rest is radiated off to infinity. The resulting exterior 
spacetime agains settle down to Schwarzschild with the Yang-Mills field in its vacuum state. However, the gauge 
potential in the slightly super-critical Type II collapse approaches its other vacuum value: w — > +1. Similarly, 
further into the Type II super-critical regime (where we should emphasize the discretely self-similar nature of the 
critical solution is no longer clear and where the black holes which form all have finite mass) we find "generalized" 
super-critical Type II collapse which can be characterized by an asymptotic value of the gauge potential of w = +1. 

It is worth commenting here that our terminology has become a bit loose. We have been speaking of Type I and 
Type II collapse well into the super-critical regime when in fact these terms originally referred to types of critical 
behavior in a very small neighborhood of the exact critical solutions. For instance, as one moves away from the Type II 
threshold, the discrete self-similarity of the exactly critical solution is no longer clearly exhibited by configurations 
which collapse to a black hole. In addition, a significant fraction of the mass present will end up in the black hole. 
This is to be contrasted with the tiny black holes that can be formed in near-critical Type II collapse. It should also 
be made clear that Type I and Type II do not refer to different types of black holes. The end-state exterior geometry 
in all of these evolutions is a Schwarzschild black hole. However, we are, in a sense, generalizing the concepts of Type I 
and Type II collapse to distinguish two kinds of collapse dynamics in the super-critical regime. 

It quickly becomes clear that these regions are well defined in parameter space and it becomes natural to try and 
find the boundary separating them. Indeed, the problem reduces to tuning a single parameter across the threshold 
between these two types of collapse dynamics. What we find is that at the threshold between these two types of 
dynamics a static solution emerges as the intermediate attractor for this new near-critical collapse. The Type II 
side of this threshold exhibits the following near-critical evolution. As the field collapses, a finite mass black hole 
is formed. The remaining field exterior to the black hole then approaches a static, non-trivial configuration. After 
some time the field then departs from the static configuration and disperses to infinity leaving the original black hole 
virtually unchanged. The gauge potential in this case approaches the vacuum value of w — + 1 and thus this side of 
the threshold is part of the generalized Type II collapse dynamics. The Type I side of this threshold has, of course, a 
very similar development initially. Again, a black hole of finite mass is formed exactly as before and is followed by the 
approach of the exterior Yang-Mills field to the static solution. However, as the field leaves the static solution, most 
of the matter collapses into the black hole with only a small portion radiating off to infinity. The gauge potential, in 
this case, approaches the vacuum state of w = — 1. Figures [l] and || show the development of these two cases overlayed 
with the static, intermediate attractor which is the exactly critical solution. 

Since we are in the super-critical regime, these new, exactly critical solutions are themselves black holes with Yang- 
Mills hair. Indeed, these attractors are the colored black holes discovered by Bizon and independently by Volkov and 
Galtsov [ p~3|Jl4| shortly after the discovery of the Bartnik-Mckinnon solutions. In general, these colored, static black 
hole solutions are described by their horizon radius (or, relatedly, their mass) Th and the number n of zero crossings 
of the Yang-Mills potential w. The solutions which serve here as intermediate attractors in the dynamical collapse 
are the lowest lying colored black holes with n = 1. Subsequent to their discovery, it was shown that these black 
holes are unstable. In particular, the colored black holes with label n have 2n unstable modes |i~8|-|2"o[] . Of these 2n 
unstable modes, n arise from spherically symmetric perturbations in the gravitational sector away from staticity and 
another n unstable modes arise from perturbations in the sphaleron sector, i.e. away from the magnetic ansatz for 
the Yang-Mills fields. Since we are solving the nonlinear evolution equations while remaining within the magnetic 
ansatz we would expect our data to excite all the unstable modes within the gravitational sector while exciting none 
of the unstable modes in the sphaleron sector. As a result, we would expect the n = 1 family of solutions to have 
a single unstable mode, thus making it relevant to critical collapse. The picture of this new critical solution in the 
black hole regime having a single unstable mode thus accords perfectly with previous understanding of both Type I 
and Type II collapse. 

This new critical behavior within the black hole regime is, in many ways, qualitatively similar to that of the 
Type I collapse. Solutions in one of these new interpolating families asymptote to one of the static, colored black 
hole solutions, Fi(r;r^). Initial data for collapsing configurations close to one of the critical solutions will approach 
Yi(r; r/i), and remain in its vicinity for some amount of time T, as measured, for example, by an observer at infinity. 
The evolution will then peel off the static solution with the remaining field either dispersing to infinity or collapsing 
and adding a finite amount of mass to the already present black hole. One can quantify the amount of time spent near 
the static solution in the same way as with the Type I collapse. The time, as measured by an asymptotic observer, 
is T s» — Aln|p — p*\, where the coefficient A is the characteristic time scale for the collapse of the unstable solution 
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Y\(r\ rh) or, in other words, the inverse Lyapounov exponent of the single unstable mode. 

With our evolutions, we confirm this relation and can calculate the value of A in the manner described in 0. For 
Figures |l| and [|, where the static solution has a black hole radius of rh ~ 0.55, we confirm the linear relation and find 
A « 4.88 as shown in Figure [| 
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FIG. 1. Time sequences of the Yang-Mills potential w for both sub- and super-critical evolutions, overlayed with a static, 
colored black hole solution. Both evolutions are calculated to machine precision, i.e. \p — p*\/p* < 10~ 15 . The solid line is 
the static solution, the dashed line is the generalized Type I solution, where part of the Yang-Mills field of the colored black 
hole collapses to a second, larger black hole, and the dotted line is the generalized Type II solution, in which the remaining 
Yang-Mills field surrounding the colored black hole disperses to infinity. The calculations shown here use the kink data defined 
in Eq.(p^). We fix the center of the pulse to be x = 2.16, and perform our parameter search on s (= p), the width of the 
pulse. These plots use the original radial coordinate, r. The initial colored black hole that forms has a radius of r ~ 0.55 on 
a domain from r = to r = 40; a fixed uniform radial grid with 3201 points is used for all calculations. In the generalized 
Type I case, the second and larger black hole forms with a radius of r ~ 1.69. Note that after the black hole forms, we do 
not plot points at locations within the horizon radius. This illustrates that our horizon excision continues the evolution only 
outside the black hole. As can be seen, the agreement between the static solution and the intermediate attractor to which the 
dynamical solutions temporarily evolve is excellent. The static solution is obtained by solving the static equations (ODEs) for 
a colored black hole of radius r — 0.55. 

Once we have established the existence of these new, colored critical solutions, it then becomes straightforward to 
map out parameter space for our choice of initial data. Figure does this for the main region of interest. 
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Figure [f| should be compared to Figure 1 in §. There, it was suggested that there was a coexistence between black 
holes exhibiting both types of critical behavior. However, that suggestion was based on evolutions which did not 
actually continue to the future of the resulting black hole. Thus, up to this point, only the region marked "dispersion" 
had been well explored. Evolutions had been performed in the black hole regime, but detailed understanding of 
what happened in collapse was largely limited to times prior to the formation of the black hole. Now, with our use 
of apparent horizon boundary conditions, we can explore more completely this super-critical region. Thus, what is 
actually present is two well-demarcated regions in the black hole regime separated by a critical line, as described 
above. 
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FIG. 2. Time sequence of the Yang-Mills potential as shown in Figure — here the waveforms are plotted using a logarithmic 
radial coordinate. 
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FIG. 3. Plot of the elapsed time (as measured by an observer at infinity) in generalized Type II collapse before the zero of 
the w field crosses r = 35. As the exactly critical solution is approached in parameter space, i.e. p — > p* , the dynamical fields 
spend more time on the critical solution; specifically, the lifetime of the static configuration is T ~ —A In \p — p*\. Since the early 
part of the dynamics is unchanged as p — > p* , the total time before the lingering field escapes and passes a finite radial value 
depends only on \p — p* |. Thus we can estimate A by measuring this total time as a function of In \p — p* \ . A straightforward 
least squares fit of this data yields A ~ 4.88. 
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FIG. 4. Plot of the phase space for the two parameter family of kink data given by (^). Each point represents a critical 
solution at a level of \p — p*\/p* ^ 10~ 4 and sits on one of the thresholds separating two of the three phases. The open hexagons 
represent Type I critical solutions which separate configurations that disperse from those that form a finite size black hole. 
We confirm that the intermediate attractor in this case is the n = 1 Bartnik-Mckinnon solution, as found in Q]. The open 
triangles represent Type II critical solutions which separate configurations that disperse from those that form an infinitesimal 
black hole. The filled hexagons represent the new black hole critical solutions separating the generalized Type I and Type II 
collapse dynamics as described in the text. These critical solutions are the n — 1 colored black holes. They are parameterized 
by their horizon radius, such that as the "triple point" is approached along the critical line, the size of the black holes goes to 
zero. 
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V. DISCUSSION 



We have presented evidence for the existence of a new type of critical phenomena within the black hole regime of 
the spherically symmetric Einstein- Yang-Mills model. Previous work on critical behavior in this model has established 
that both Type I and Type II critical solutions emerge in sufficiently general families of initial data. However, that 
earlier work was unable to accurately evolve to the future of black hole formation. Using a different coordinate system 
and horizon excision techniques, we are able to evolve well into the black hole regime. We find that these two different 
types of critical collapse can be generalized to types of collapsing spacetimes with a distinguishing characteristic being 
the asymptotic vacuum value of the Yang-Mills potential. These two "phases" of black hole formation are separated 
by a critical line in phase space such that dynamically forming black hole solutions exactly at criticality are the static, 
n = 1 colored black holes jl3|,[l4]]. This, of course, fits nicely within our conceptual framework that critical solutions 
are intermediate attractors of co-dimension one. 

It is worth emphasizing a few observations at this point. First, it is important to note that the mass of the final 
Schwarzschild black hole in the evolution will be discontinuous across this new critical line in Figure ^. (In fact, indi- 
cations of jumps in the black hole mass-spectrum in the super-critical regime were present in earlier calculations plf .) 
Naively, this might be somewhat surprising since it means that nearby configurations in the space of initial data can 
have, in principle, very different sized black holes as their end-states. 

Further, we note that when the colored black holes were originally investigated, it was found that the nth static 
colored black hole is parameterized by its horizon radius, r^. In addition, as r/j — > 0, the black hole solution approaches 
the corresponding nth Bartnik-Mckinnon solution. In the dynamical context discussed here for the n = 1 solutions, 
one notices an interesting manifestation of this result. The colored black hole attractors separating the Type I and 
Type II regions in Figure ^[ are likewise parameterized by their horizon radius in such a way that as the "triple 
point" is approached, ry t — > 0, and the Bartnik-Mckinnon solution emerges as the relevant attractor separating Type I 
collapse from dispersion. 

Finally, it is natural to ask if other collapsing systems might exhibit similar behavior. Indeed, one can simply look 
in the literature on static, hairy black holes and point to a number of systems that one would conjecture should 
have a "phase diagram" similar to Figure |4|. Such theories should include Einstein- Yang-Mills coupled to a dilaton 
(EYMD) and the Einstein-Skyrme system. However, there are other examples of systems which exhibit both Type I 
and Type II critical behavior. These include a massive scalar field and (it is conjectured) a massive, charged scalar 
field H . These latter systems are particularly intriguing in light of the results presented here because though they do 
exhibit Type I and Type II critical behavior, no-hair theorems suggest that these systems do not have static black hole 
solutions. It should thus be of interest to investigate the nature of the super-critical regime for these, and possibly 
other, systems. 
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APPENDIX: 



We collect here some of the details involved in deriving the equations of motion used in the paper, and to present 
a somewhat more general framework for future work on the Einstein- Yang-Mills system (EYM). The action for our 
model is 



S = 



R 



16ttG g 2 



_pa pa^iv 



(Al) 



where F^ v is the Yang-Mills field strength tensor. As usual, Greek indices range over the four spacetime dimensions 
and Latin indices indicate group indices. Spacetime indices are raised and lowered by the metric g^ v while we will 
not lower group indices. Thus, repeated upper group (Latin) indices will be summed over. 

The equations of motion are found by varying the action with respect to the fields. Varying with respect to the 
metric yields the Einstein equations 
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16ttG G ^ ~ T ^ 

= ^{2F; x F\^- 1 -g, v F^F^). 

Varying with respect to the connection A® gives 



(A2) 



(A3) 



With the general equations in hand we are in a position to make some simplifying assumptions. To begin, we assume 
that we have an SU (2) gauge field. In addition we will assume spherical symmetry. The most general parameterization 
of the spacctime metric is then 



ds 2 



(-a 2 + a 2 /3 2 )dt 2 + 2a 2 (3dtdr + a 2 dr 2 + b 2 r 2 dQ 2 



(A4) 



where the metric coefficients a, f3, a and b will depend on the temporal and radial coordinates t and r and dil 2 is the 
usual metric on the unit two-sphere. Likewise, the most general parameterization of the gauge connection is now 

A = ur r dt + vr r dr + (wr e + wr' f ')d8 + (cot 8 t t + wt^ - w t ) sin 8d(j). 

where the coefficients u, v, w and w will all depend on t and r. The r a are the spherical projection of the Pauli spin 
matrices and form an anti-Hermitian basis for the group SU(2) satisfying [r a , r b ] = e ahc T c with a, b, c g {r, 8, </>}. We 
note that this connection is invariant under a gauge transformation of the form U = e^(*' r ) TT '. The field strength 
derived from this connection is 



F = r r {v-u')dt Adr 

+ [(w - uw)dt + (w' - vw)dr] A (r 6 'd8 + sinddcj)) 
+ [(w + uw)dt + (w' + vw)dr] A {r^dd 
-(1 - w 2 - w 2 )r r d8 A sin <9d0 



r sm 



(A5) 



With these assumptions, the equations of motion for the Yang-Mills fields can now be written in first order (in 
time) form as 



n 

P 



f3U+ -$ 
a - 

Hp + -Q 

a 



+ uP-v[pP+-Q)+ -p^wil - w 2 - w 2 ) 
V a > b z r z 



- uU + 



v ( pn + - 

V a 



b 2 r' 



-w(l — w 2 — w 2 ) 



(A6) 
(A7) 



where we have used the definitions 



Y' = wU - wP 

OL 

Y=-(w$- wQ) + (3(wU - wP) 



II = — \w — uw — 3(w' — vw)] 
a 

$ = w — vw 

P = — \w + uw — f3(w / + vw)] 
a 

„7.' 



Q = w + vw 

b 2 T 2 



Y = 



2aa 

We also have evolution equations for $ and Q: 



(v - u') 



-n + /3$ +uQ-v -P + /3Q 



_ 2aa^ r 
— W tttt;Y 



Q= -P + /3Q 
.a 



u$ + v 



-n + /?$ 

a 



+ w 



b 2 r 2 
2aa, 



b 2 r 2 ' 



(A8) 
(A9) 

(A10) 

(AH) 
(A12) 
(A13) 
(A14) 



(A15) 
(A16) 
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which follow from diff erentiatio n of (All) and (A.13) with respect to time, and elimination of explicit time derivatives 
using the definitions (A10-A14). 

Finally, the relevant Einstein equations are the evolution equations for the metric components and the components 
of the extrinsic curvature, together with the Hamiltonian and momentum constraints, as follows: 



a = —aaK r 
b = -abK e 



{aP)' 

P < hY 
~ ( rb > 



(A17) 
(A18) 



k r r = (3K r ' r + aK r r K--(-) - —^ 

a \ a I arb 



K% = (3K 8 ' g + aK%K 



a 



1 



(rb)' 



(rb) 2 a(rb) 2 \ 



( arb 



A-nGa [S-p- 2S r r ] 
(rb)' ) + AnGa [S - p - 2S%] 



(A19) 
(A20) 



arb (a)^r£>^(a^^) J 



AK r r K% + 2K e l = WnGp 



(A21) 



rb 



K\) - K% = AnGj 



(A22) 



Here, p is the energy density, j r is the momentum density and S r r , S g and S^^, are the stress energy components 
projected onto our spacelike hypersurface. Explicitly, we have 



V {b 4 r 



1 f AY' 2 (1 - w 2 - w 2 ) 2 



b 4 r 4 



b 2 r 2 a? 



[q 2 + $ 2 + p 2 + n 2 ] 



(A23) 



b r - V 



AY 2 

" & 4 r 4 



(1 - w 2 



b 4 r 4 



b 2 r 2 a 2 



[q 2 + $ 2 + p 2 + n 2 ] 



(A24) 



AY 2 
Ag 2 \¥r~ 4 



{l-w 2 



(A25) 



(A26) 



g 2 ab 2 r' 



(n$ + pq) . 



(A27) 



As can be seen, we have chosen to write things in a fairly general form. Indeed, this is the most general EYM 
theory we could consider in spherical symmetry. We can thus consider more general matter configurations, together 
with a greater variety of coordinate systems than have been considered to date. Subsets of these equations have been 
investigated before in different contexts. For example, static versions of these equations have been studied in order 
to understand particle-like solutions and colored black holes. In addition, evolved a version of these equations in 
polar (K e g = 0), areal (b = 1) coordinates with the additional assumption on the Yang-Mills field that the connection 
A^ a was purely magnetic. 

For our work here, we will also make this "magnetic ansatz." More specifically, we will assume that the electric 
charge density, Y, is identically zero. In addition to this ansatz, we also make the following gauge choice: v = 0. 
Making this gauge choice leaves the connection invariant under a gauge transformation of the form U = e^^ T . From 
the definition for Y in equation (A14) we see that the function u, by these assumptions, depends only on t. However, 
the residual gauge freedom implies that u is arbitrary up to a function of t. We can thus completely fix the remaining 
gauge freedom by choosing u — as well. A consequence of these choices is that the remaining coefficients of the 
gauge connection will no longer be independent and one of them (we choose w) can be set to zero without loss of 
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generality. Thus, in the context of this particular problem, the magnetic ansatz has effectively set all the fields but w 
to zero. 

As discussed in the body of the paper, we want to be able to evolve the system for long times to the future of any 
black hole formation. To this end, we choose a coordinate system with maximal slicing: K — K r r + K 6 g + K^^ = 0. 
We will retain the choice of areal, or radial, coordinates (b = 1) so that the coordinate r is immediately related to the 
area of origin-centered two-spheres. 

Given all these assumptions, the equations which we must solve simplify considerably. The evolution equations are 



IT 

w 



a 

a 



+ —w(l - w 2 ) 



ft 



and the constraint equations are 

w' = $ 



a' 2 
a r 
\~a 2 



2a 



1 



2ra' 



AirGa {S - 3p) 



K e 9 ' 



2r 

3 -k, 

r 



2 

4ttG 



ra 3 K e e2 + AirGrap 



g 2 ar 2 



where the matter stress-energy terms are given by 

a 2 {l-w 2 ) 2 



S-3p 



2g 2 r i 
a 2 (l-w 2 ) 2 



g 2 r 2 



1 



Ag 2 r 4 



2g 2 r 2 



$ 2 +W 



($ 2 + W 



Finally, we have an algebraic condition for the non- vanishing component of the shift vector 

(3 = arKg 9 , 



(A28) 
(A29) 
(A30) 

(A31) 
(A32) 

(A33) 
(A34) 

(A35) 
(A36) 

(A37) 



which follows from (A18) with b(r,t) = 1. 
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